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Because scaling symmetries of the Euler-Lagrange equations are generally not variational symme- 
I , tries of the action, they do not lead to conservation laws. Instead, an extension of Noether's theorem 

^ ■ reduces the equations of motion to evolutionary laws that prove useful, even if the transformations 

' are not generalized symmetries of the equations of motion. In the case of scaling, symmetry leads 

to a scaling evolutionary law, a first-order equation in terms of scale invariants, linearly relating 
kinematic and dynamic degrees of freedom. 

QThis scaling evolutionary law appears in dynamical and in static systems. Applied to dynami- 
cal central-force systems, the scaling evolutionary equation leads to generalized virial laws, which 
. linearly connect the kinetic and potential energies. Applied to barotropic hydrostatic spheres, the 

scaling evolutionary equation linearly connects the gravitational and internal energy densities. This 
i— 1 1 implies well-known properties of polytropes, describing degenerate stars and chemically homoge- 

i-C . neous non-degenerate stellar cores. 

1 PACS numbers: 45.20.Jj, 45.50.-j, 47.10.A-, 47.10.ab, 47.10.Df, 95.30.Lz, 97.10.Cv 
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I. SCALING SYMMETRY NOT GENERALLY A SYMMETRY OF THE ACTION 

Action principles dominate physical theories because they admit transformations among dynamical vari- 
ables and exhibit common structural analogies across different systems. If these transformations are symme- 
tries of the action, then by Noether's theorem, they give rise to conservations laws that reduce the number 
of degrees of freedom. This relationship between symmetries of the action (variational symmetries) and con- 
servation laws is central to Lagrangian dynamics. But, even if these transformations are not symmetries of 
the action, they nonetheless lead to useful Noether's identities. Although equations of motion do not require 
Lagrangian expression, we apply this identity to transformations that are not variational symmetries; in 
particular, to scaling symmetry, which is generally a generalized symmetry only of the equations of motion 
(Section II). Variational symmetries and generalized symmetries both reduce the equations of motion to first 
order, but in different ways: 

• Variational symmetries imply conservation laws, first integrals of the equations of motion. 

• Scaling symmetry generally implies only an evolutionary equation, which reduces the equations of motion 
to first order in scaling invariants. 

Applied to dynamical systems of bodies interacting via inverse power-law potentials, these scaling evolu- 
tionary equations are generalized virial theorems (Section III) . Applied to self-gravitating barotropic spheres 
in hydrostatic equilibrium (Section IV), the scaling evolutionary equation will be an analogous first-order 
equation between scaling invariants. In this way, scaling evolutionary equations will illuminate the physical 
consequences of scaling symmetry. 

So as to focus on scaling evolutionary equations, we relegate our Lagrangian formulation of barotropic 
hydrostatics to Appendix A, and needed stellar thermohydrodynamics to Appendix B. A following paper [l| 
will derive the well-known properties of polytropes and of homogeneous stellar cores from Section IV of the 
present paper. 

We do not consider applications to quantum field theories @ , involving the symmetry of the vacuum as 
well as the Lagrangian, which lead to important quantum anomalies and to topological symmetries generated 
by topological charges. 
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II. NOETHER'S THEOREM EXTENDED TO NONINVARIANT TRANSFORMATIONS 



A. Noether's Identity Implies Either Conservation Laws or Evolutionary Equations 

We start with a system of particles described by the Lagrangian C(t, qi, q\) and action S — j C(t, qi,qi)dt, 
where the dot designates the partial derivative d/dt with respect to the independent variable and the Einstein 
summation convention is assumed. Under an infinitesimal point transformation 5(t,qi),Sqj(t,qi) generated 
by 5t ■ d/dt + Sqi ■ d/dqi, velocities and Lagrangian transform locally as 

d5qi dSt 

SC = [St ■ d/dt + 5 qi ■ d/d qi + Sq t • d/dqi]C = C5t + (dC/dq t )5q % + (d£/dqi) - q~\ , (1) 

where d/dt = d/dt + q\ ■ d/dqi + 9i • d/dqi is the total derivative. The canonical momentum and energy 

p l {t,q l ,q l ) := dC/dqi , E(t,qi,qi) := q\{dC/dqi) - C (2) 

have total derivatives 

d ^ = dL/dq i -V i , -^ = C + V t -q\, (3) 
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in terms of the Euler-Lagrange variational derivative T>i := dC/dqi — d(dC/d(ji)/dt. Since 



4;(Pi s <li) = {dC/dqi - Vi)5qi + p { ■ d ^ 
at at 



dESt) ,. Nr d(St) 



the Noether charge 



has total derivative 



G := £ ■ 5t + pi ■ (5qi - q\St) = -E5t + piSqi 
d{8t) 



dt 



Vi ■ (Sqi - qidt) — SC -T>i- (Sqi - qtSt) , 



(4) 



(5) 



(0) 



where SC := SC + C ■ (dSt/dt) is the change in Lagrangian at a fixed point. 
The variation in action between fixed end points is 



5S 12 



dtdC 



dt 



^_ £ .^ +v ..( 6q ._ 4 . St) ] =G (i)_ G (2)+ f dt 
at dt i ,h> 



r ^ c fdh dC\ 



(7) 

after integrating the term in d{8t)/dt by parts. The action principle asserts that this variation vanishes for 
independent variations 5qt,8t that vanish at the end points. It implies the Euler-Lagrange equations and 
dh/dt — —dC/dt. On-shell, where the equations of motion = hold, 



SSl-2 = 



SC dt = G(l) - G(2) 



dt 



(8) 
(9) 



This is Noether' s Equation, an evolutionary equation for any generator in terms of the Lagrangian trans- 
formation that it generates. It expresses the equations of motion as the time derivative of the Noether 
charge. 



B. Variational Symmetries Imply Conservation Laws 



The most general and important applications of Nocther's Equation are to variational symmetries and 
to generalized symmetries of the equations of motion, which both preserve the stationary action principle 
8S\2 = but reduce the equations of motion to first order in different ways. 

Variational symmetries preserve the action 5S\2 = because 8C = or dB/dt, the total derivative of some 
gauge term B(t, q). Noether's Equation (d/dt)(G — B) = —T>i ■ (Sqi — qiSt) conserves G — B on-shell, where the 
equations of motion hold. This original Noether's Theorem, identifying conservation laws with space-time 
variational symmetries, is familiar in the symmetry of central-force systems (|13[) under time translations and 
spatial rotations, leading to conservation of energy E and angular momentum I: 

E := (f 2 + r 2 9 2 )/2 + V(r) , l:=mr 2 . (10) 

An illustration of how gauge terms B(t : c\i) conserve G — B, rather than the Noether charge G, is the 
many-body system of particles with interparticle forces that depend only on the relative separations — 
and relative velocities — q^ . This system admits the infinitesimal boost transformations 

8<ii = 5v-t , St = , 8V = , 8K = SC = P ■ Sv , (11) 

where M, P, K are the total mass, momentum, and kinetic energy. Because boosts change all the momenta, 
the charge G = (P • v)i is not conserved. Instead, Noether's equation gives the conservation law (P — 
A/R) ■ Sv = 0, or MR = P, for arbitrary infinitesimal Sv. While boosts change the total momentum P, 
the center-of-mass moves with velocity R. This familiar center-of-mass theorem follows directly from boost 
symmetry, irrespective of the internal forces. It is paradigmatic for distinguishing between the effects of 
internal and external forces in many-body system. 

Conversely conservation laws imply invariance of the Lagrangian modulo a possible gauge term, so that 
variational symmetries imply conservation laws. 
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C. Generalized Symmetries Reduce the Equations of Motion to First Order 



The remainder of this paper will deal with generalized symmetries of the equations of motion, which are 
not variational symmetries and do not lead to conservation laws, but to transformation equations. These 
generalized transformations, nevertheless, reduce the equations of motion to first order. For example, a 
nonrelativistic particle, with radial, transverse and angular momenta 

p,. = (r-p)r/r 2 , p f = p - p r = -r x l/r 2 , 1 := r x p , (12) 

and energy 

E(r,p) =K-V(r) := p 2 /2m - V(r) , p 2 = p 2 r + l 2 /r 2 , (13) 

in an inverse-power central potential V(r) ~ l/r". The Lagrangian for this system £(r, p) = K + V(r) := 
p 2 Jim + V(r) admits the scaling symmetry 

5t = 0t , Sr = r , Sp = (1 - (3)p , SK = 2(1 - (3)K , 5V = -nV, (14) 

generated by 

G nA := -0Et + r • p, (15) 

provided j3 = 1 + n/2, so that the Lagrangian is homogeneous of degree —n^JjC = 2(1 — fj)C — —n£. Since 
5(rP /t) = 0, all distances scale with time as r, ~ t 1 ^ , as shown in Table I Q. 
The scaling evolutionary equation 

dG nA /dt = SC=(1- n/2)C , (16) 

linearly connects the kinetic and potential energies to the time derivative of the single particle virial A := 
r • p = rp r . This scaling evolutionary equation is the virial equation 

d(rp r )/dt=(l-n/2)(K + V) + (l-n/2)(K-V) = 2K + nV . (17) 

Only for zero-energy orbits E ~ K + V ^ in &n inverse-cube force n = 2, is the virial A conserved and 
scaling a variational symmetry. 

The first-order orbit equation for the scaling invariant 

rdO/dr =pt/p r = l/^r 2 p 2 /l 2 - 1 = 1/ y/2mr 2 [E - V(r)]/l 2 - 1, (18) 
can be solved by quadratures 



9(r)=0 Q + dr/{r^2mr 2 [E - V{r)]/l 2 - l}. (19) 

J ro 

(In the Kepler Problem V(r) = —k/r, the integrals reduce to elementary functions and the bound orbits are 
ellipses 

r{6) = all - e 2 )/[l - esin(0 - 6 % (20) 



TABLE I: Period- Amplitude Relations and Virial Theorems for Inverse Power-Law Potentials V ~ l/r 



n 


System 


Period-amplitude relation t ~ r 1+n ' 2 


Virial theorem 


-2 


isotropic harmonic oscillator 


period independent of amplitude 


{K) = 


(V) 


-1 


uniform gravitational field 


falling from rest, e.g., z = gt 2 /2 


(K} = 


(V)/2 





free particles 


constant velocity r ~ t 


(K) = 





1 


Newtonian potential 


Kepler's Third Law t 2 ~ r 3 


{K) = 


-(V)/2 


2 


inverse-cube force 


t~r 2 


{K) = 


-(V) 



5 



of eccentricity e := y/T— I 2 /mka and semimajor axis a := —k/2E.) From this first-order equation and 
angular momentum conservation, the temporal evolution is 



dt = (mr 2 /l)d8 = mdr/^2m(E - V)(r)) - (l/r) 2 , (21) 

so that 



t(r)=t + dr/\/2r 2 [E -V(r)} - (l/m) 2 (22) 
completes the integration of the central-force problem. 

D. Even Transformations That Are Not Symmetries Lead to Useful Evolutionary Equations 

For the same system, the radial translation 

6t = f3t , 5r = r/r , 5p = -p t /r , (23) 

generated by 

G nB :=p-Sr=p r (24) 
is not a generalized symmetry. The radial translation evolutionary equation 

dp r /dt = l 2 /mr 3 - dV/dr (25) 

is the radial equation of motion. 

Generalized symmetry or not, the first-order evolutionary equations (17) for rp r and (25) for p r are useful 
expressions for the radial equation of motion. While the virial A and radial momentum p r are not conserved, 
their time averages (dA/dt) = (dp r /dt) = in any bounded crgodic system, so that, for time averages 

scaling transformation imply 2(K) = — n(V), the usual virial law generalized to arbitrary inverse-power 
potentials; 

radial displacement transformation imply (Z 2 /m)(l/r 3 ) = —n(V/r), I / 0, useful for relativistic cor- 
rections to noncircular Coulombic (n = 1) orbits [J]. 

The next section considers energy-conserving many-body systems C=C(ri,Ti), for which the scaling evo- 
lutionary equation is a generalized virial law. Section IV, treats the hydrostatics of barotropic spheres, for 
which the Lagrangian C(r, H, H') depends explicitly on the radial variable r. Instead of a first integral, In 
both these examples, instead of a conservation law, scaling symmetry implies an evolutionary equation, a 
first-order differential equation between scaling invariants [5j, |6| , linearly relating the the "kinetic" term K 
to the "potential" term V in the Lagrangian C = K — V . 

III. SCALING EVOLUTIONARY EQUATIONS 

A. Mechanical Evolutionary Equation is a Generalized Virial Law 

The previous derivation generalizes to a nonrelativistic many-body system with particle coordinates r^. 
The scale transformation 

St = /3-t , 6r l = n , 8{d/dt) = -p ■ (d/dt) , 6Ti = (l-0).Ti , SK = 2(1 - p)K (26) 
is generated by the Noether charge 

G n :=-(3-Et + A , (27) 
where A :— J^p; • r^ is the many-body virial. 



G 



If the pairwise potential energies are inverse powers Vij ~ r,^™ of the interparticle distances r.y := |rj — ry), 
the interparticle potentials are homogeneous in their coordinates r(dVij/dr) = —nVij, SVij — — nV«. Scaling 
is still a generalized symmetry of this many-body system, provided n = 2(j3 — 1), f3 = 1 + n/2. Then 
SV = — (1 — /3)r • VV, 6 A = (1 — 0)A and the total energy and Lagrangian are homogeneous functions of 
their arguments, scalar densities of weight —n, 

5E = -nE , 5C = -nC , 5C = 6C + f3£ = (1 - n/2)C. (28) 

Because energy is conserved, Noether's equation (9) implies the scaling evolutionary equation 

dG 

—^ = -(l + n/2)£ + A = S£ , A = (1 + n/2)E+ (1 -n/2)£ = 2K + nV . (29) 
at 

This is a generalized virial law linearly relating the time derivative of the virial A to a linear combination of 
the nonrelativistic kinetic energy K and the power-law potential V . 

For periodic or long-time averages in bounded ergodic systems, (A) = 0, so that the virial theorem relates 
time averages 2(K) = —n(V) . Table I tabulates these period-amplitude relations and generalized virial 
theorems for orbits in the five important inverse-power-law potentials n = —2, —1, 0, 1, 2. Only for inverse 
cube forces V(r) ~ 1/r 2 would dynamical scaling reduce to a symmetry of the action, and the Noether 
charge Gi = — 2{K + V)t + A be conserved. For potentials more singular than 1/r 2 , there are no bound 
states. 



B. Scaling Evolutionary Equation in Classical Electrodynamics 

Noether's equation applies to continuous Lagrangian systems (fields) as well as discrete systems. In this 
case, r,t are independent variables. If f , G = Ex B/47rc, T, U are respectively the electromagnetic force 
density, momentum density, momentum flux tensor, and energy density, then momentum balance reads 

dG/dt + V -T + i = . (30) 

From this follows, the evolutionary equation : 

d(r-G)/dt + V-(T-r)-U + r-f = , (31) 

an electromagnetic analogue of the mechanical virial law (31). When time-averaged, this becomes an elec- 
tromagnetic virial theorem 



C. Scaling Evolutionary Equation in Classical Conformal Field Theory 

In any relativistic field theory, space-time scaling (dilatation) symmetry leads to the familiar evolutionary 
equation 

% , (32) 



dxf* 



where d/dx^ is the four-dimensional divergence, G^ is the dilatation current, and 6£J is the trace of the 
energy- momentum tensor @, 0j@l- If this trace vanishes, the dilatation charge is conserved, implying con- 
formal symmetry. 

The most familiar example of conformal symmetry is Laplace's equation in n spatial dimensions. In two 
dimensions, conformal symmetry implies the Cauchy-Riemann equations, so that any analytic function is a 
solution of Laplace's equation. In higher dimensions, conformal symmetry implies the conservation laws asso- 
ciated with translations, rotations, dilatations, and spatial inversions. Ignoring charges, the electromagnetic 
field is conformally invariant. 

These familiar examples from conservative systems recall how Noether's equation leads to useful evolu- 
tionary equations, whether or not scaling symmetry is broken. The chief purpose of this paper remains to 
consider the hydrostatic equilibrium of barotropic spheres, which differ by being static, not Hamiltonian 
(since the radial coordinate r is not ignorable), and derive from a Least Energy variational principle, instead 
of a Least Action principle. 
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IV. SCALE-INVARIANT BAROTROPIC STARS 
A. Mechanical Structure of Barotropic Stars 

The hydrostatic structure of barotropes depends only on mass continuity and pressure equilibrium, 

dm/dr = Airr 2 p , -dP/dr = Gpm/r 2 , (33) 

or the second-order equation 

1 d /r 2 dP\ 

■\^-(-?f)=-4nGp , (r 2 H'y+4irGr 2 p(H)=0 , (34) 
H dr \ p dr / 

where ' := dj dr. In terms of the entropy H (r) = J dP/ p and specific gravitational force dV/dr = g := 
Gm/r 2 , the equation of hydrostatic equilibrium reads 

d{H + V)/dr = , (35) 

so that is Poisson's equation. We consider only barotropic stars for which the local equation of state 
is P = P(p). and the density p{r), specific internal energy E(r), specific enthalpy H{r) = E + P/p, and 
thermal gradient V(r) := dlogT / dlog P are implicit functions of the gravitational potential V(r). 
This static structural equation is the Euler-Lagrange equation of the Lagrangian 

C(r,H,H') = 4irr 2 [-{H') 2 /%it + P{H)], , (36) 

derived from a minimal energy variational principle in Appendix A. It describes the radial evolution of a 
static barotropic sphere. Because the radial coordinate r appears in C(r, H, H'), the energy per radial shell 
E(r, H, H') = — 4irr 2 [(H') 2 /8itG + P{H)] increases while moving outwards in radius. 



B. In a Simple Ideal Gas, Scale Invariance Requires a Constant Entropy Gradient 

Polytropes are barotropic spheres in which the polytropic index n := dlog p/ dlog H or polytropic exponent 
1+1/n := d log P/d log p are constant (Appendix B). The pressure, specific energy, specific enthalpy, enthalpy 
gradient and central pressure at any point are 

P/P c = (p/p c ) 1+1/n , E = n(P/p) , H = (n + l)(P/p) , dlog H/dlogP = l/(n + 1) . (37) 

The hydrostatic structure is polytropic in zero-temperature (degenerate) stars and nearly constant in 
convective stars and in stars starting out on the hydrogen-burning, zero-age Main Sequence (ZAMS), where 
the chemical composition and energy generation are homogeneous: 

White dwarfs and neutron stars: nonrelativistic and extreme relativistic degenerate stars; polytropes 
of index n=3/2 and 3, respectively. 

ZAMS stars in convective equilibrium: with vanishing gravithermal specific heat C* and uniform en- 
tropy density. These are n=3/2 polytropes. 

ZAMS stars in radiative equilibrium: At zero age, our Sun was a chemically homogeneous star of 
mean molecular weight p, — 0.61, well-approximated by the Eddington standard model n=3 polytrope 
throughout its radiative zone, containing 99.4% of its mass. 

Because energy generation was extended but not uniform, our ZAMS Sun would be better approximated 
globally by a slightly less standard n=2.796 polytrope 0. Still better nonpolytropic fits would obtain 
by including both nonuniform energy transport and corrections to Kramers opacity: radiative transport 
in the pp-burning lower main sequence 0.11 < M/M© < 1.2 gives the exponent £ = 0.57 in the M — R 
relation R ~ M*; convective transport in the CNO-burning upper main sequence 2 < M/Mq < 20 gives the 
exponenet £ = 0.8 US El-) 

Because our present Sun is chemically inhomogeneous and has convective zones, it is far from being 
polytropic; its best polytropic fit, with index n = 3.26, is poor Q. 



C. Scaling Symmetry Implies a First-Order Equation in Scaling Invariants 

Following Chandrasekhar [ljj], we define homology variables 

u := d\ogm/d\ogr = 3p/p , v n :— — d\og (P/p)/d\ogr , w n := — dlogp/dlogr = n(r) ■ v n , (38) 

where p = Sm/Airr 3 is the average mass density interior to radius r and n(r) :— d\ogp/d\og(P/p). The 
central boundary condition is 

u(0) = 3 , v(0) = , (dv/du) = -5/3n . (39) 

The mass continuity and hydrostatic equilibrium equations (|33j) become 

d\ogu/d\ogr = 3 — u — n{r)v , d\ogv/d\ogr = u — 1 + v — d\og [1 + n(r)]/d\ogr , (40) 

and will be autonomous only when the index n(r) is constant. 

In polytropes, the constant index n and gradient V = l/(n + 1) makes both these equations 

du/d log r = u(3 — u — nv n ) , dv n /d\ogr = v n (u — 1 + v n ) (41) 

autonomous, so that they can be written as the characteristic differential equations 

d log v n d log u 



u — 1 + v n 3 — u — nv r 



= d\ogr (42) 



for the homology invariants u, v n jl7| . 
The infinitesimal scale transformation 



Sr = r , 5H = -dJ n H , 8H' = -(1 + L) n )H' (43) 
is generated by the Noether charge 

G n :=-E-r+p-5H = -r :i [^ + AitP{H)]-r 2 ^ r -uj n H . (44) 

The structure is scale-invariant if and only if P — Kp 1Jrl / n , so that H ~ P/p, P ~ H n + l and 

^ \ 

7 )' 

n — 1/ 

(45) 

making the gravitational potential and specific internal energy both homogeneous of degree — 2Cj n and 

5C = -2w n C , 5E = -2w n E . (46) 

The radial derivative 

dG n /dr = 5C + C + V r -u n - d{Hr) jdr , (47) 
obeys the scaling evolutionary equation 

dG n /dr = <f„£ , (48) 

connecting the gravitational and internal energy densities, just as the virial law connected the potential and 
kinetic energies. 

For such polytropes, we introduce dimensionless units 

i-=r/a , 9 n :=H/H c = (p/p c )^ n , (49) 

and the dimensional constant 

a 2 := { -^^Kp 1 J n - 1 = (n+l)/4nG-(P c /p 2 c ) , (50) 



9 



where p c is the central density and 

P c l Pc :=Kp 1 J n , H c :=(n+l)P c /p c . 
The included mass, mass density, average mass density, and gravitational acceleration are 

m = 4ir Pc a 3 • (-£X) . P = Pc-K , p = Pc-(-K/0 > 5 = 47rp c a 2 (-^) 



(51) 



(52) 



Poisson's equation (|34|) . combining mass continuity and hydrostatic equilibrium, takes the dimensionless 
Lane-Emden form 



d fgj d °n 



Suppressing the subscript n on 8 n and 0' n := d9 n /d^, the homology variables are [l^ 



(53) 



u := dhgm/dlagr = 3p(r)/p = -£6 n /6' 
u/v n =9 n+1 /e' 2 

Extracting the dimensional constant C 



v n :=-d\og(P/p)/dlogr = -Z0'/9 , p/ Pc = -W'/e , 



UV n = £ 



Inn— I 



\n-l 



, (uv n ) 



l/(n-l) _ <rw„+l/_fl/ 



(54) 



-H^/G, the Lagrangian, Hamiltonian, and Noether charge are 



c/c = e 2 



q/2 



gn+1 



-2w„rV"^ri 



(™) 



2 n+ V 
1/2 



1 2 n + 1 



gn+1 



n + 1 



) + U)r, 



2 n + 1 

-2oj" 



n+1 



U 2 n+1 J+^nU V„ 



For n = 5, scaling symmetry is a variational symmetry and the Noether charge 



Ok-?- 



5^,3/2 



(™ 5 ) 3/2 



+ (m, 5 3 )V2 



is conserved. Otherwise, the Noether charge evolves according to ((48 



n + 1 



(55) 



(56) 



(57) 



which reduces to the Lane-Emden equation (|53[) . This evolutionary equation describes the growing ratio 
between local internal and (negative) gravitational energy densities 



gn+l 



/(n + 1) 



#' 2 /2 



(58) 



as the local energy density changes from entirely internal at the center, to entirely gravitational at the stellar 
surface. 



V. CONCLUSIONS 



We have extended Noether's Equation connecting variational symmetries to conservation laws to general- 
ized transformations of the Euler-Lagrange equations. The resultant evolutionary equations are useful, even 
when the transformations are not generalized symmetries. But, when they are generalized symmetries, they 
reduce the Euler-Lagrange equations to first-order equations between invariants. 

For scaling symmetries, the evolutionary equation takes a special form connecting linearly the kinematic 
and dynamic parts of the Lagrangian. For nonrelativistic systems with inverse-power law potentials, the 
scaling evolutionary equation is a generalized virial law, linearly relating the kinetic and potential energies. 
For hydrostatic systems obeying barotropic equations of state, the scaling analogous evolutionary equation 
is a linear relation between the local gravitational and internal energies. From this scaling evolutionary 
equation, in the following article [l|, we derive all the properties of polytropes. 
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Appendix A: LAGRANGIAN FORMULATION OF BAROTROPIC HYDROSTATICS 

Stellar structure generally depends on coupled equations for pressure equilibrium and heat transport. 
Only if the heat transport leads to a local barotropic relation P = P(p) can the hydrostatic equations be 
considered independently, without reference to the thermal structure. 

1. Mass Continuity and Hydrostatic Equilibrium 

In a self-gravitating isolated system in local thermodynamic equilibrium, a barotrope held at zero external 
pressure has thermodynamic potential energy, the work needed to adiabatically extract unit mass, or specific 
enthalpy H{p) = E + P/p. Barotropic energy conservation, dH := dP/p, makes the specific enthalpy a 
more natural state variable than the specific internal energy E, pressure P, or density p. The equation of 
hydrostatic equilibrium —dP/dr = Gmp/r 2 := pg is then 

-dH/dr = dV/dr = g = Gm(r)/r 2 , (Al) 

describing how this local specific enthalpy or extraction energy H(r) depends on the local gravitational 
potential V(r). Integrating, we have the energy conservation equation 

H(r) + V(r) = -GM/R , r <R , (A2) 

where the zeros of the gravitational potential and specific enthalpy have been chosen at infinity and at the 
spherical surface, respectively. 

Because the gravitational potential obeys Poisson's equation 



±±( r i dV 
r 2 dr V 

the specific enthalpy obeys the second-order equation 



V 2 H + AirGp(H) = . (A4) 



To implement the equation of hydrostatic equilibrium, we need a local barotropic relation P(p), or P(H), 
p(H), which is determined by the thermal stratification of the static matter distribution in local thermody- 
namic equilibrium, and by a central boundary (regularity) condition (dP/dr) = = (dp/dr) . Near the 
origin, 



p*p c (l-Br 2 ) , m (r)*^f-(l-^Br 2 ) 



p{r) := 3m(r)/4^r 3 » p c (l - ^Br 2 ) w p 2 J 5 p 3 / 5 . (A5) 

In terms of the homology variables w := — rflogp/rflogr, u := dlogm/dlogr for the mass density and 
included mass, dw/du w —5/3. 

2. A Constrained Minimum Energy Principle for Hydrostatic Equilibrium 

In a static, self-gravitating sphere of mass M and radius R, the Gibbs free energy 

W:=E-TS + PV = n + U , (A6) 
in terms of the gravitational and internal energies 

i-M j-R 

= - / (Gm/r)dm , U = - PdV , (A7) 
Jo Jo 
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where p, E, and —Gm(r)/r are the mass density, specific internal energy, and gravitational potential, re- 
spectively. In the Eulerian description, the radial coordinate is r, the enclosed volume is V = 47rr 3 /3, and 
the enclosed mass m(r) is constrained by mass continuity dm{r) = pdV. The Gibbs free energy 

rR rR. 

W= C(r,m,m')dr = - Awr 2 [Gmp/r + P{p)]dr , ' := d/dr , (A8) 
Jo Jo 

is the available work when expanding the sphere adiabatically at fixed external pressure. The Lagrangian C 
is the Gibbs free energy per radial shell dr. 

The constrained minimum energy variational principle (TlLIT^ for hydrostatic equilibrium is that the Gibbs 
free energy be stationary (5W = 0) under adiabatic deformations in specific volume SV P — d(Aitr 2 5r) / dm 
that vanish on the boundaries and satisfy the mass continuity constraint m! = Aitr 2 p. This minimum energy 
principle has the equation of hydrostatic equilibrium 

V r := Gm/r 2 + H' = Gm/r 2 + dP/pdr = , (A9) 

as its Euler-Lagrange equation, with mass continuity as a constraint. This equation is scale invariant if the 
specific enthalpy H scales as ml . 



3. An Unconstrained Variational Principle 

Using Poisson's equation to incorporate the mass continuity constraint, the gravitational energy is 

n = - (V' 2 /2)ATrr 2 dr , (A10) 



so that the second-order Lagrangian (used in Section IV) 

C(r, H, H') = Anr 2 [-H' 2 /8ttG + P(H)} (All) 

is unconstrained and has Euler-Lagrange equation (IA4[) . The canonical momentum and Hamiltonian are 

p := dC/dH' = -r 2 H'/G = -m , H(r, H,p) = -Gp 2 /2r 2 - A-Kr 2 P{H) , (A12) 

and the canonical equations are 

dH/dp = H' = -Gp/r 2 , dU/DH = -p = m = Anr 2 p . (A13) 

Spherical geometry makes the system non- autonomous, so that dH/dr — —dC/dr — —2C/r and m'(r) vanish 
only at large r, when the geometry approaches being planar. 



Appendix B: STELLAR THERMOHYDRO DYNAMICS 



The structure of luminous stars depends upon the coupling between hydrostatic and thermal structures 
through an equation of state P — P(p, T, which generally depends on the local temperature and chemical 
composition. But, ignoring evolution, the matter entropy is locally conserved, so that steady-state stars 
are in local thermodynamic equilibrium (LTE). In a fluid held in pressure equilibrium at constant external 
temperature, the specific Gibbs free energy H — TS = — V(r) is a minimum. In hydrostatic equilibrium, the 
density p(r), specific internal energy E(r), specific enthalpy H(r) = E + P/p, specific entropy and thermal 
gradient V(r) := <i log T/<i log P depend implicitly on the gravitational potential V(r). 

In the Second Law of Thermodynamics 

TdS = dQ = dE + Pd{l/p) , (Bl) 

E, p can be written as functions of temperature and pressure. Clever use of thermodynamic identities then 
leads to [13, [Hi 

TdS=C*dT , dS = c P (V - V ad )d log P , (B2) 

where the gravithermal specific heat C* := dS/d\ogT = P(l — V a( j/V) depends on the specific heat cp and 
the adiabatic gradient V a d := (<91ogT/<91ogP)s. This expression relates the local thermal gradient V(r) to 
the local entropy gradient dS(r)/dlog P — cp(V — V a d), which generally contains both gas and radiation 
entropies and varies in stars in radiative equilibrium. 
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In a simple ideal gas , the equation of state, specific internal energy, specific enthalpy and adiabatic 
exponent are 

P/p=—T , E = c v T , H = c P T , (d log P/d log p) s = c P /c v := 7 » (B3) 
M 

where 7?. is the universal gas constant, /x is the molecular weight, and cv, cp — cy + are the 
specific heats at constant volume and at constant pressure. From the Second Law of Thermodynamics 
(B2), the specific entropy and thermal gradient of a simple ideal gas are 

dS = c v d log P-cpd log p , S = c v log (P/p 7 ) , V = d log P/d log P = (7- l)/ 7 . (B4) 

In an ideal gas supported by both gas pressure P gas = ~R.pT/ fi and radiation pressure P ra d = aT 4 /3 
, the total pressure is 

P = Pgas + Pad = /3P+ (1 - P)P = (KpT/p)[l + (1 - P)/p\, (B5) 

where the pressure ratio 

Pad/Pgas := (1 -P)IP= MR) ■ S rad , (B6) 

is proportional to the specific radiation entropy SVad = 4aT 3 /3p. The total specific entropy 

S = S g + S iad = (TZ/p)- [log (T 3 / 2 /p) + 4(1 -13)/ (3} . (B7) 

is therefore constant for a cool monatomic ideal gas (/3 = 1) or for a radiation-dominated (supermassive) 
star(/3«0). 

Bound in a polytrope of index n , an ideal gas has constant thermal gradient, gravithermal specific 
heat, and entropy-pressure gradient: 

V = 1/0+1) , c* = c P (l - V ad /V) ,dS/d\ogP = c P (V -V ad ) <0. (B8) 

Except for this, a polytrope's thermal structure is unconstrained and still depends on the unspecified 
heat transport mechanism. 

According to Schwarzschild's minimal entropy production criterion, convective stability requires that a star's 
specific entropy stay constant in convective equilibrium and increase radially outwards in radiative equilib- 
rium. This makes barotropic stars of mass M extremal in two respects: the central pressure is minimal in 
for a given radius P; the central pressure and temperature are maximal for a given central density. Because 
stellar evolution is driven by developments in the core, these bounds drive stars toward uniform entropy in 
late stages of evolution 14 1 . 



Acknowledgments 

SAB acknowledges support from the Millennium Center for Supernova Science through grant P06-045-F 
funded by Programa Bicentenario de Ciencia y Tecnologia de CONICYT and Programa Iniciativa Cientifica 
Milenio de MIDEPLAN. 



[1] S. Bludman and D. C. Kennedy (????), arXiv:1008.4281. 

[2] M. E. Peskin and D. V. Schroeder, An Introduction to Quantum Field Theory (Perseus Books, 1995), section 
19.5. 

[3] L. D. Landau and E. M. Lifshitz, Mechanics (Addison- Wesley, 1976), 3rd ed., section 10. 

[4] J. Schwinger, J. L. L. Deraad, K. A. Milton, W. Tsai, and J. Norton, Classical Electrodynamics (Perseus Books, 
1998), sections 3.3, 8.5. 

[5] P. J. Olver, Applications of Lie Groups to Differential Equations (Springer- Verlag, 1993), 2nd ed., ISBN 3-540- 
94007-3, exercise 5.35. 



13 



[6] W. Blumen and S. Kumei, Symmetries and Differential Equations (Springer- Verlag, 1989). 
[7] J. C. G. Callan, S. Coleman, and R. Jackiw, Annals of Physics 59, 42 (1970). 
[8] S. Coleman, in Aspects of Symmetry (Cambridge University Press, 1985), chap. Chapter 3. 
[9] S. A. Bludman and D. C. Kennedy, Astroph. J. 525, 1024 (1999), figures 2, 3; Table 1. 
[10] R. Kippenhahn and A. Weigert, Stellar Structure And Evolution (Springer- Verlag, 1990), ISBN 3-540-50211-4, 
figure 22.2, Table 20.1. 

[11] C. J. Hansen and S. D. Kawaler, Stellar Interiors: Physical Principles, Structure, and Evolution (Springer- Verlag, 

1994), section 1.2; Figures 7.4, 7.5. 
[12] S. Chandrasekhar, An Introduction To The Study Of Stellar Structure (University of Chicago, 1939), chapters 

III, IV. 

[13] H.-Y. Chui, Stellar Physics (Blaisdell Publishing Company, 1968), section 2.12. 
[14] A. Kovetz, Mon. Not. R. Astr. Soc. 144, 459 (1969). 

[15] W. E. Boyce and R. C. DiPrima, Elementary Differential Equations and Boundary Value Problems (John Wiley 
and Sons, 2001), seventh ed. 

[16] D. W. Jordon and P. Smith, Nonlinear Ordinary Differential Equations (Oxford University Press, 1999), 3rd ed., 
problem 2.13. 

[17] This evolutionary equation is important elsewhere in mathematical physics. In population dynamics, with logr 
replaced by time t, it becomes the Lotka-Volterra equation for predator/prey evolution [TBI \W^. The uv n cross- 
terms lead to growth of the predator v n at the expense of the prey u, so that a population that is exclusively 
prey initially (v n = 0,u = 3) is ultimately devoured u 0. For the weakest predator/prey interaction (n = 5), 
the predator takes an infinite time to reach only the finite value v n = 1. For stronger predator/prey interaction 
(n < 5), the predator grows infinitely v n ~ ¥ oo in a finite time. 



